function rest=myfun1(t,y,a,b,m)

rest=a*m-a*y;